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In this paper we study the finite size scaling for low energy 
excitations of S = 1 and S — 2 Heisenberg chains, using the 
density matrix renormalization group technique. A crossover 
from 1/L behavior (with L as the chain length) for medium 
chain length to 1/L 2 scaling for long chain length is found 
for excitations in the continuum band as the length of the 
open chain increases. Topological spin S — 1/2 excitations 
are shown to give rise to the two lowest energy states for both 
open and periodic S = 1 chains. In periodic chains these two 
excitations are "confined" next to each other, while for open 
chains they are two free edge 1/2 spins. The finite size scaling 
of the two lowest energy excitations of open S — 2 chains is 
determined by coupling the two free edge S = 1 spins. The 
gap and correlation length for S = 2 open Heisenberg chains 
are shown to be 0.082 (in units of the exchange J) and 47, 
respectively. 

PACS Numbers: 75.10.-b, 75.10.Jm, 75.40.Mg 



Since Haldanc pointed out the difference between in- 
teger and half integer Heisenberg spin chains Q, the 
richness of the physics in the Heisenberg model has at- 
tracted intensive attention. It was conjectured |jj that 
the ground state of integer spin chains has an exponen- 
tially decaying correlation and a nonzero energy gap for 
excitations, in contrast to the ground state for half inte- 
ger spin chains with a power law decay correlation and 
a zero energy gap. Like valence bond solid picture 0, 
free edge spins are expected to exist in both integer and 
half integer open Heisenberg spin chains with site spin 
5 > 1/2 [|L For 5 = 1 spin chains, the field theory 
prediction [0] for the continuum band and magnetic sus- 
ceptibility of S = 1 chains agrees with experiment ||; 
the energy gap is found both numerically |^| and experi- 
mentally . The numerical calculation on effects of free 
edge spins for S > 1/2 chains ^ also agrees with pre- 
vious analytic studies on edge states dti. For 5 = 1 
chains, the finite size scaling for energy is expected to 
be 1/L 2 (L is the chain length) for the multi-magnon 
states with periodic boundary condition (PBC) ||, and 
the same 1 / L 2 scaling is also anticipated for open bound- 
ary condition (OBC) chains with one or several magnons 
||, while the first excitation energy scales exponentially 
for PBC chains ||. In contrast, different scaling for- 
mula have been used in studying 5 = 2 chains |^o| [L4| . 
Therefore, a systematic study of finite size scaling for the 
energy spectrum of integer spin chains at various chain 
length is highly needed to obtain a precise gap value for 
5 = 2 spin chains. 



In this paper, we study numerically the finite size scal- 
ing for low excitation energies in the continuum band 
for 5=1 and 5 = 2 Heisenberg open chains, which 
exhibits a i behavior for medium chain length and 
dependence for long chain length. The crossover be- 
tween medium length behavior and long length behav- 
ior is demonstrated by numerical results for 5=1 open 
chains. Such a crossover reconciles the numerical scal- 
ing found recently Jll|,|l2| for 5 = 2 spin chains and the 
theoretical scaling for chain length L — > oo. The finite 
size scaling of energies for the two lowest states of 5 = 1 
Heisenberg chains is given and topological spin 1/2 exci- 
tations are shown to be responsible for them. We show 
explicitly that the two spin 1/2 edge states are "confined" 
to each other in PBC chains. The finite size scaling of 
the two lowest excitation energies in 5 = 2 OBC chains 
is also determined by their free edge 5=1 spins. The 
gap and correlation length for 5 = 2 Heisenberg chains 
are shown to be 0.082(3) (in units of the exchange J) 
and 47(3), respectively. We also point out that a reason- 
able scaling for the first excited state (the bottom of the 
continuum band) of 5 = 2 PBC chains is of exponential 
type in analogy with 5=1 PBC chains. 

We apply the density matrix renormalization group 
(DMRG) method |15j to calculate the lowest energies for 
Heisenberg chains using the Hamiltonian: 
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(The above equation is for OBC. The Hamiltonian for 
PBC chains is obtained by adding the term Sl • Si to 
Eq.(l).) We use the DMRG infinite chain algorithm to 
calculate the low energy spectrum for both OBC chains 
and PBC chains fl5|| . By keeping 300 optimized states, 
we calculate the energies Eg 1 for the ground state and 
.E-f 1 for the first excited state in 5 = 1 PBC chains up 
to length L = 36. The largest truncation errors in cal- 
culations are 9 x 10~ 8 for Eq 1 and 5 x 10~ 7 for -Ef 1 , 
respectively. By keeping 200 optimized states, we calcu- 
late the energies Eq 1 for the ground state and E® 1 for 
the first excited state of 5 = 1 OBC chains up to length 
L = 42. The largest truncation errors in the calculations 
are 2 x lO" 11 for Eg 1 and 1 x 10~ 13 for E? 1 , By keep- 
ing 150 optimized states, we calculate also the energies 
E§ 1 , E® 1 , and E^ 1 for the low-lying states of total spin 
Sl ot = 2, 3, and 4, respectively for 5 = 1 OBC chains 
up to length L = 100. E^ 1 and E® 1 are then calculated 
for chain lengths up to L = 250. The largest trunca- 
tion errors occurred in these calculations are 1 x 10 -12 , 
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2 x 1(T 10 , 1 x 1(T 9 , and 6 x 1(T 9 for Eg 1 , Eg 1 , Eg 1 , 
and Eg 1 , respectively. By keeping 250 optimized states, 
we calculate the energies Eg 2 , Eg 2 , Eg 2 , and Eg 2 for 
5 = 2 OBC chains up to length L = 150, and they are 
corresponding to the lowest energy states with total spin 
5* ot = 0, 1, 2, and 3, respectively. The largest trunca- 
tion errors in the calculations are 1 x 10~ 6 , 2 x 10~ 7 , 
4 x 10- 8 , and 7 x 10~ 8 for Eg 2 , Eg 2 , Eg 2 , and Eg 2 , 
respectively. The analysis of these numerical results is 
straightforward, and we will show them in detail in the 
following. 

We use the following relations for the two lowest ener- 
gies in S — 1 Heisenberg chains (to be explained later): 

Eg^eoL + a^e-^/VL, 

E^^eoL + A.+a^e-^^/VL, f . 

Eg 1 = e (L - 1) + A 6 + a™ e -C*-i)/fc [Z) 
Eg 1 = e (L - 1) + A 6 - ia^e-^" 1 )/^. 

for PBC (the first two equations) and OBC (the last two 
relations), respectively. It is more convenient to intro- 
duce 

h{L) = -(EP- - e Q L)VL = -a^e-Wi, 
h(L) = (EC 1 - e L - AJy/L = of 1 e- L /^, 
h{L-l) =-(£ 01 - eo (L-l)-A b ) = -oOi e -(i-i)/€x, 
U{L-l) = Eg'-eoiL-D-A, = -ffie-W**. 

(3) 
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FIG. 1. The logarithm of the energy expressions in Eq.(3) 
for the two lowest energy states of S = 1 spin chains as func- 
tions of chain length L. 

Our numerical results are plotted in Fig.l as bxfi(L) 
v.s. L in comparison with Eq.(3). From fitting we obtain 
the site energy eo = —1.401484(1), correlation length 
£i = 6.00(4), energy gap Ai = 0.4104892(2), and edge 
state energy for OBC chains A b = -0.1931661(1). The 
constants are of similar precision for £i, a^ 1 = —1.33(1), 
afVao 1 = -0-33(2), and ag 1 = -0.50(2). 

The interpretation of OBC results is obvious: Eg 1 
and Eg 1 are the energies of the singlet and triplet com- 
binations of the two spin 1/2 edge states, respectively. 



Hence the ratio of the constants is — g and the cou- 
pling constant J eff ~ er^- 1 )/* §,|Q. It is very in- 
teresting to note that a^ 1 /ag 1 ~ —1/3 in Eq.(2), which 
means that the singlet and triplet combinations of the 
two topological spin 1/2 excitations are also responsible 
for the two lowest states for PBC chains, with a cou- 
pling constant J e ff ~ e~ L ^/s/L. An important ques- 
tion is where are these topological excitations located. 
For OBC chains they are exactly pinned at the edges. 
One might think that these two states are pushed away 
as far as possible from each other for PBC chains. How- 
ever, this is not true: they are "confined" , sitting next 
to each other |]l7| . These two excitations will interact 
with each other in "both directions" , a direct coupling 
of the order of J which is responsible for the Haldane 
gap, and an indirect coupling via virtual spin-wave ex- 
citations around the spin chain which is proportional to 
J eff ~ exp(-L/£)/VX @. 

The exponential term e~ L ^l\f~L in Eq 1 and Ef 1 is 
in agreement with previous exact diagonalization results 
Q. By studying chain length longer than the one used 
in exact diagonalization, we have shown that Eq 1 and 
E[ 1 have the same scaling behavior as for the spin-spin 
correlation P,pl|, instead of a net exponential e^ L ^. We 
point out that Ey 1 is the bottom of the one magnon band 
and has a zero kinetic energy since w(k) ~ (k — 7r) 2 fflJR, 
therefore its scaling is different from the 1/L 2 behavior 
expected for other single or multi-magnon states in PBC 
chains as well as the magnon states in OBC chains. 
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FIG. 2. The logarithm of the energy spacings (Eq.(4)) 
for the three lowest energy states in S = 2 spin chains as 
functions of chain length L. Diamonds are for \n(Eg 2 -Eg 2 ), 
and crosses are for \n(Eg 2 — Eg 2 ). 

Now we turn to S = 2 chains. The asymptotic energy 
differences for the three lowest states in S = 2 OBC 
chains are: 

Eg 2 -Eg 2 = ag 2 e-^- x y^, 
Eg 2 -Eg 2 = 2a° 2 e-( i - 1 )/« 2 , (4) 

where the correlation length £2 = 47(3) and constant 
ag 2 = 0.060(5). The DMRG data are plotted in compar- 
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ison with the above equation in Fig. 2 as \n(E^ 2 — Eq 2 ) 
and ln(£'^' 2 — E® 2 ) v.s. L — 1. These asymptotic expres- 
sions come from the coupling of two spin 5 = 1 edge spins 
in 5 = 2 OBC chains ||,|^]. From these energies we have 
calculated the site energy for 5 = 2 Heisenberg chains 
e = -4.76125(5) by E$ 2 (L + 2) - E$ 2 (L) analysis jjj). 
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FIG. 3. The excitation energies for the bottoms of one, 
two, and three magnon bands of medium length 5 = 1 OBC 

-(n-l)Ai](i-l) v.s. 
20 to 100. 



chains (Eq.(5)). We plot [E% 1 - E[ 
l/ln(L — 1) for medium lengths L 
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FIG. 4. The excitation energies for the bottom of one 
magnon band of medium length and long length S = 1 OBC 
chains. We plot (E§ 1 - E? 1 - Ai)(L - 1) v.s. 1/ ln(L - 1) for 
medium lengths L — 20 to 100, and for long lengths L — 100 
to 250. The crossover from (E§ x - E? 1 - Ai)(L - 1) » 
" l '•" ' 1)) for medium lengths (dashed line) to 



a? 1 



(l + 6?7ln(£ 



Ai)(L 



74.7/(L - l) 2 for long lengths 



(E^ - E\ 

(full line) M is evident. The bottom of the magnon dispersion 
E? 1 - E? 1 = y'A 2 + (2.49tt/L) 2 § is also plotted (dotted 
line). 



As for excitations in the continuum, we first describe 
the medium length scaling behavior for integer spin OBC 
chains. For S = 1 chains, we use the following energy 
expression for medium length: 



E 



oi 



E" 1 = (n - l)A a 



a: 



ox 



b oi 



1 (1 ~' ln(L-l) 



(5) 



where n — 2,3,4 correspond to one, two, and three 
magnon band bottoms, respectively. We plot in Fig. 3 
(L-lXE^-J^-Ai), ( J L-l)( J Bf 1 -£; l 01 -2A 1 ), and 
(L-1)(E2 1 -E^ )1 -3A 1 ) v.s. l/ln(L-l) for S = 1 OBC 
chains with medium chain length, L — 20 to 100, which 
is more than ten times the correlation length £i. The 
ratios of the constants are a^/a^ 1 ~ 4, a ox ja® 1 » 9, 
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b^ 1 . For Jong chains, L > 100, the ex- 
pected scaling term 1/L 2 



and b 01 



_ for E° x - ^f 1 shows up. 
In Fig.4, using again the plot (L - ^(Ef? 1 - E? 1 - A x ) 
v.s. l/ln(i — 1) for chain length L = 20 to 250, we 
show the crossover from the behavior given by Eq.(5) for 
medium chain length to E® 1 - E® 1 = A t + J^lyi @ 
for long chain length. We point out that the crossover 
occurs within a rather narrow range in 1/ In L plot in 
Fig.4. For S = 2 OBC chains, we use medium chain 
length expressions similar to the ones for 5=1 chains, 
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1 v ln(i - 1 
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where the gap turns out to be A2 = 0.082(3) and con- 
stants a° 2 = 2.8(2) and bf 2 = 5.9(4). We plot in Fig.5 
(L - 1){E§ 2 - E$ 2 - A 2 ) v.s. l/ln(i - 1) for chain 
length up to L = 150. As we show in Fig.5, the medium 
chain length behavior extends to chain length L = 150 
and bigger. For the chain lengths we have studied, which 
is more than three times the correlation length, the lead- 
ing term 1/L in Eq.(6) is still dominant. In analogy to 
5=1 chains HJ^l , the long length asymptotic expression 



for S = 2 OBC chains is expected to be -jj as well [[12| 
and the crossover between the medium length and long 
length behavior should occur for longer chain lengths. 

One way to explain the 1/L scaling in the medium 
range is to consider the energy dispersion A(fc) = 
v/A 2 + (vk) 2 for a magnon with wave vector k. This ex- 
pression has also been suggested for such gapped systems 
using nonlinear sigma models p(J . For OBC chains with 
big length L, one can expand the (vk) 2 term and get the 
1/L 2 scaling. So it is in agreement with our calculation 
and previous studies 0. The leading finite size scaling 
term would be 1/L only when A 2 <C (2ttv/L) 2 . The 
shape of the curve a/A 2 + (2irv/L) 2 v.s. L is consistent 
with Eg 1 - Lf 1 (see Fig.4), but it can not fit Eg 1 - £f 1 
for the medium size well, for instance, L ~ 20 to 90 for 
5 = 1 OBC chains. We suggest that further studies on 
the medium size behavior of 5 = 1 and 5 = 2 chains are 
needed, especially because it is relevant to experiment on 
chains with typical doping (of the order of one percent). 
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FIG. 5. The excitation energies for the bottom of one 
magnon band of medium length S = 2 OBC chains (Eq.(6)). 
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We plot (E' 3 
L = 50 to 150. 
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-A 2 )(X-1) v.s. l/ln(L-l) for lengths 



The leading order scaling 1/L for medium chain length 
has been found numerically in recent studies pr| , ^2[ for 
5 = 2 Heisenberg chains under various boundary con- 
ditions. Given the decreasing of the gap equal to, or 
more rapid than, the inverse of chain length 1/L, our 
calculation gives a lower bound 0.082(3) for the gap of 
S = 2 chain. In a previous work ]12j, a similar lower 
bound 0.085(5) was also given for the gap using a differ- 
ent boundary condition [^l|,^2| . Since the lower bound of 
the gap is finite, the existence of a nonzero gap in S = 2 
Heisenberg chain has been proved by numerical calcula- 
tions, and therefore, Haldane's conjecture is verified also 
for S = 2 chain. Moreover, the calculated gap is also in 
agreement with an earlier calculation in Ref. fll4| and a 
recent experiment p3j. We emphasize that the gap cal- 
culation in the present paper disagrees with some other 
studies Jl^,|ll| . This discrepancy may arise from the dif- 
ferent scalings used, such as 1/L 2 scaling [lj| for the gap 
of short S = 2 PBC chain, where an exponential type 
scaling is more appropriate. 

To conclude we have performed a systematic study of 
finite size scaling behavior for the energy spectrum of in- 
teger spin chains (S — 1, 2), using the DMRG technique. 
This analysis shows that the lowest excitations are due 
to the "edge" states for both OBC and PBC chains. As 
for excitations in the continuum a crossover from 1 / L be- 
havior for medium chain length to 1 / L 2 behavior for long 
chain length has been demonstrated. Using this analysis, 
a reliable estimate for the gap and correlation length in 
S = 2 chains has been obtained. 
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Theoretical Physics, Trieste. 
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